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§ . Abstract 
O 

p • We consider a micromaser model to study the influence of Dicke superradiance in the 

Ch ' context of the one-atom maser. The model involves a microwave cavity into which two- 



level Rydberg atoms are pumped in pairs. We consider a random pump mechanism which 
allows the presence of at most one pair of atoms in the cavity at any time. We analyze the 
' differences between the present system, called the Dicke micromaser, and an equivalently 

pumped conventional one-atom micromaser. These differences are attributed to the Dicke 
cooperativity in the two-atom system. We also show that the two-atom Dicke micromaser 
is equivalent to a one-atom cascade two-photon micromaser. With the introduction of a 
one-photon detuning, the present theory further describes a true two-photon micromaser. 
We discuss in detail the role of one-photon detuning in the mechanism of a one-atom 
two-photon micromaser. This leads us to point out that the two-atom cavity dynamics 
can be verified by a proper scaling of the results from an equivalent one-atom two-photon 
micromaser. 
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1. Introduction 

The first experimental observation of maser action witli just one atom in tlie so-called 
micromaser [1] renewed interest in the subject extensively as it opened possibilities of 
observing cavity-QED effects which are purely of quantum mechanical origin. Indeed, 
one of the results of the much discussed Jaynes-Cummings interactions [2], the quantum 
revival of the oscillations in the atomic population, has been observed in the micromaser 
device [3] . The micromaser device [1] consists of a high-Q superconducting cavity cooled 
down to sub-Kelvin temperatures into which Rydberg atoms in the upper of the two 
masing levels are pumped randomly but one at a time. A clever velocity selector decides 
a fixed flight time r through the cavity for every atom. The repetition time tn, r and tcav, 
the duration in which the cavity is empty of atom, satisfy in = T+tcav with r <^ tcav This 
is the basis of the one-atom maser theories [4-6] with the further assumption of strictly 
single-atom events. 

However, the above arrangements in the experiments [1,3] cannot eliminate overlap 
of flight times through the cavity of successive atoms completely, which of course, is less 
than 1%. This has been one of the obstacles in realizing certain predictions of one-atom 
maser theories [4,5] such as number states of the cavity radiation field. Another limiting 
factor is the dissipative mechanism acompanying the coherent atom- field interaction [6]. 
The results in [6] show that the effects of the dissipative mechanism can be minimized by 
increasing the Q factor of the cavity and by decreasing its temperature. But the pump 
mechanism has not been able to eliminate traces of two atom events in one-atom maser 
dynamics. Hence it has become necessary to explain the effect of such two-atom events 
on the photon statistics of the cavity field. Various models incorporating multiple atomic 
events have been proposed [7-13]. The basic difference between such models and one- 
atom maser theory [4-6] is the Dicke atomic cooperative effect [14,15] acompanying the 
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multiple atomic events. Thus, in order to understand the limitations of maser theories 
based on strictly single atom events [4-6], it is necessary to quantify the effect of atomic 
cooperation. For this purpose we propose the following model. 

We consider a fixed number of atoms entering the cavity with a constant flight time 
through it, but with random time gaps between successive such events. The cooperative 
nature [12] of the interaction of the atoms with the cavity field shall be fully refiected in 
the photon statistics. In this paper we consider atoms being pumped in pairs into the 
cavity. We limit the number of atoms to two just for the sake of simplicity which can 
be easily generalized to a larger number. Both the atoms are taken to be in their upper 
states when they enter the cavity. The arrival of the pair at the cavity is assumed to 
be Poissonian. We thus have tn — t + tcav with Ir — 1/R. We assume r to be a fixed 
duration for every pair of atoms, and hence, tcav is random. R is now a Poisson average 
of R pairs of atoms pumped into the cavity in one second. Due to cooperative interaction 
of the pair of atoms with the cavity field, we will find its dynamics with R pairs of atoms 
pumped into the cavity per second in the present case being different from the one-atom 
maser pumped with 2R atoms into the cavity per second. This is due to the superradiance 
effects [14,15] which will be manifested in the steady-state photon statistics. For these 
reasons, we call the present system a "Dicke micromaser" . 

It is rather a difficult task to achieve such a model experimentally. However, we will 
show in the following that the present model is closely related to a one-atom two-photon 
micromaser [16-18] . This is a key result of this paper. Two-photon micromaser action with 
one atom has been experimentally demonstrated [16]. Thus the model we are considering 
to show the atomic cooperative effects is not beyond experimental verification. 

The organization of the paper is as follows. In the section 2, we present the microscopic 
Hamiltonian and the Dicke states involved in the micromaser action. The equations of 
motion are derived in section 3. We derive and discuss the photon statistics in section 
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4. The present system is compared with the two-photon micromasers in section 5. We 
conclude the paper in section 6. 



2. The Model 

A pair of Rydberg atoms in the upper of their two levels a and b enter the cavity at 
t — 0. We assume that the atomic transition frequency is in resonance with the cavity 
mode frequency cu. Both the atoms take the time r to travel across the cavity. During 
this time, the evolution of the system comprising of the two atoms and the single eigen 
mode of the microwave cavity is governed by the Hamiltonian 



as we consider resonance condition between atomic transition frequencies and the cavity 
frequency. Hi is the Jaynes-Cummings Hamiltonian [2]. 

Hj = gj:{Sra + S-a^) (Ic) 

i=l 

We assume that the atom-field couphng constant g is the same for each atom. The 
atomic operators obey the relations [S^, S^] = 2Sf6ij and a{a'^) is the photon annihilation 
(creation) operator obeying the commutation relation [a, a^] = 1. In the frame rotating 
at Hq, we have H — Hj. 

It is convenient to represent the two-atom system by a Dicke system as both the atoms 
interact with the same cavity field simultaneously. In doing so, we however, neglect the 
dipole-dipole interaction as the microwave cavity dimensions are much large compared to 



(la) 



where 



Ho = u;(S^ + Si + a^a) 



(16) 
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particle wavelengths. Defining the collective operators = J2'^=iSf, Eq.(lc) takes the 
form 

Hi = g{S+a + S'a^) (2) 

which is clearly invariant by atomic permutation. Both the atoms are in their respective 
upper states at t — and hence, the state representing the atomic system \a,a) is also 
invariant under atomic permutation. So it remains in a symmetrical state at all times. 
These are known as Dicke states with maximum cooperative number J — 1 (we have a 
two-atom system) and these states are isomorphous to angular momentum states with 
principal angular momentum J = 1. These states, represented by |J, M), are generated 

by 



S-\J,M) = ^J{J + M){J - M + 1)\J, M - 1) 

and 

(3) 

S+\J,M) = ^J{J + M + 1){J - M)\J,M + 1) 

The level | J, M) is nondegenerate and has atomic energy Mu. As — J < M < J, this has 
2J +1 states, and in the present case, the number of states is three. Since J is a constant 
here, we represent |J, M) = |M). 

Thus, the pair of two-level atoms can be represented by a three-level system of the 
Dicke states \M — 1), \M — 0), and \M — —1), having energy a;, 0, and —u; respectively. 
This allows us to write the Hamiltonian in Eq.(2) in the matrix form 



Hi 



( 9\2{n-l) \ 



gp{n-l) g^2n 

V ^V2n ) 



(4) 



where n represents the photon number of the cavity field. Eq.(4) resembles the case for 
a two-photon process in a three-level system [19]. Following the method in [19], we can 
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obtain the eigenvalues and eigenstates of the interaction of the Dicke states with the cavity 
radiation field. The eigenvalues 

Ao = 0, 



A+ = g^2{2n-l\ (5) 



and its states are, respectively, 

1 0, n) = xi") 1 1 , n - 2) + 1 0, n - 1) + z^^ | - 1 , n) , 

|+, n) = 1 1, n - 2) + |0, n - 1) + 4"^ I - 1, n) , (6) 
I-, n) = |1, n - 2) + )|0, n - 1) + 4"^ | - 1, n), 
where |l,n), |0, n) and | — l,n) are the composite Dicke and field states, and 



-V2n-1' ~ hn-r ^^"^ 



(„)_ (n)__J_ in) _ / ^ .^.x 

" V4n-2' " V2' ' " V4n-2' ^^^^ 



(„)_ ,,(«)_ _L ~W _ / ^ f7.^ 

~ V4n-2' ~ ' ~ V4n-2- ^^'^ 

The states |+, n), |— , n) and |0, n) can be called as dressed states of the interaction of the 
Dicke states with the cavity field. 



3. Derivation of the equations of motion 

Using the dressed states in Eq. (6) as the basic states, the equation of motion 

p=-z[H,p] (8) 
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can be derived easily [19]. p represents the composite atom-field system having the initial 
conditions 



where P„ is the photon distribution function of the cavity radiation field. In writing 
Eq.(6), we have neglected the influences of the atomic as well as the cavity reservoirs 
on the dynamics. These effects can play a crucial role in the one-atom micromaser [6]. 
However, these influences have been shown to be negligible [6] for a cavity having very 
high Q and at very low temperatures and a low atomic pump rate. Eq.(8) represents such 
a system. 

The time-dependent solutions for the density matrix elements can be written as 



By inverting the matrix which expresses the dressed states as a linear superposition of 
the Dicke states, we can obtain the Dicke states |i, n), i = 1, 0, —1, in terms of the dressed 
states. Using these relations and Eq.(lO), we obtain the density matrix elements in the 
Dicke state basis. We follow these steps to get the P„ at t = r, which is given by the trace 



(9) 




in) in) -i{X^-^-\f)t 



(10) 



hj =0,+, -. 



+1 



P,(t)= (M,n|p|M,n). 



(11) 



M=-l 



We thus have 



n-2 



(12) 



where 





2^{n + 2){n + l) 



^(4n + 3)(4n + 6) 
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+ 



2(n+ 1) 
{An + 6) 



,(n + 2) 



)n 



(13a) 



erV) = + (4"^^^)^) - xr^^xr^^i?e.(e-(^r^-^^""^^)^), (136) 



^2n(n-l) / ) . („) 



(2n - 1) 



+^^x('^)xi")i?e.(e-(^?'-^-')-). (13c) 

The change in P„ at t = r is then 6Pn = Pn{'T) — Pn where P„ is the photon distribution 
function at the time of a pair of atoms entering the cavity. For a time At such that 
T <^ At <^ tp, where tp is the cavity photon hfetime, we have 

APn = SPnRAt (14) 

where R is the number of pairs of atoms passing through the cavity in one second 

The coarse-grained time derivative due to gain from the atomic interaction is given by 



dPn, 
dt 



\gain 



R({Qt^ - 1)P„ + Qt:^Pn-i + Qt^Pn-^. (15) 



Eq.(15) represents the dynamics during time r. During the time lapse tcavi between the 
flights of two successive pairs of atoms, the cavity field interacts with its reservoir, and 
its equation of motion is given by [6] 

Pn\ioss = 2(n + l){nth + l)nPn+i - '^n{n + nth + '^nnthjPn + 2nKnthPn-i (16) 

where k — {2tp)~^ is the cavity bandwidth, and fith is the thermal photon present in the 
cavity. Under the coarse-graining assumptions, the complete equation of motion combines 
Eqs.(15) and (16) additively. This assumption has been seen to be valid [20] for a random 

9 



input of atoms into the cavity, which is the case we are studying in this paper. We thus 
have, 



4. Steady-state photon statistics 

The equation of motion for P„, given by Eq.(17), involves Pn-2, and this makes it 
difficult to get a steady-state solution by the method followed in case of the one-atom 
micromaser [6]. However, following Risken [21], Eq.(17) can be recast in a tri-diagonal 



an expression for Tn in the form of matrix continued fractions. Eq.(17) directly gives 
an expression for Pi in terms of Pq) and Pq can be determined from the normalization 
condition. This completes the determination of P„. But numerical evaluations of P„ go 
out of control due to extremely slow convergence of the matrix continued fractions. 

Hence, wc follow the method given below, which we find computationally efficient. 
First, we set n at a value n — rimax, say, and write all the P„ in a vector 




(17) 



matrix equation involving the two-component vector 




provides 



^ — {Pq, Pi, P2, 




(18) 



This vector obviously has Umax + 1 elements. We get n. 



"max 



equations involving 



Pq, Pi, P2, 



,Pnmax froKi Eq.(17). With the normalization condition 



(19) 



n=0 



we get the equation 



MX = I 



(20) 



where M is a {rimax + 1) x {rimax + 1) matrix and / is another vector having again rimax + 1 



elements 



/ = {0,0,0, 



(21) 
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The solution 

X = M-^I (22) 

gives all Pn for < n < Umax- 

We analyze the photon statistics of the cavity field by numerically evaluating its first 
and second moments, that is, 

<n>= J2 (23) 

n=0 

proportional to the intensity of the cavity field and its normalized variance 



< > — < n >^ 



(24) 

V < n > 

i> = 1 for a coherent state field, and thus, v < 1 indicates the nonclassical nature of the 
cavity field. The process is repeated for a higher value of Umax and the corresponding 
< n > and v are compared with those obtained from the earlier value of n^nax- If they 
agree to a desired accuracy, the process is stopped, and we have the numerical values of 
the photon distribution function P„, as well as the average < n > and variance v. 

The results are displayed in Figs. 1-3. For the convenience of describing the micromaser 
action, we define the pump parameter 

D = VNqt (25) 

where N — R/2k is the number of pairs of atoms that pass through the cavity in a photon 
lifetime {2k.)~^. The patterns in the variation of < n > and v as D{(x r) changes, can 
be understood from the dynamics [22] of a collection of two-level atoms interacting with 
the radiation field of a single mode cavity with each atom being coupled to the cavity 
field by the Jaynes-Cummings interaction [2]. For a two-atom case with initial condition 
J = 1, it has been shown in [22] that the dynamics is controlled by a d.c. term with 
a prefactor of the order of Pn/n and terms oscillating at the first and second harmonic 



of the Rabi frequancy 2gsjn + ?>/2 of the interaction Hamiltonian in Eq.(lc) and having 

11 



prefactors proportional to P„ and Pn/n respectively. In the single-atom case [23] (Jaynes- 
Cummings model), the well known dynamics is controlled by only one term having the 
prefactor P„ and oscillating at the Rabi frequency 2g\/n + 1 of the Jaynes-Cummings 
interaction [2]. These differences make the present results different, in general, from the 
one-atom micromaser action. Figs. 1-3 show the differences clearly where we compare the 
results for the present Dicke micromaser for N — 100 with the one-atom micromaser 
action [4-6] with a pump rate N = 200. The numerical values of the two pump rates 
make the total number of atoms that pass through the two cavities having the same Q 
equal, and thus, justifies the comparison between the two systems. Fig. 2a shows that the 
superradiancc nature of the interaction in Eq.2 makes the threshold values of D lower 
compared to the one-atom micromaser. The differences between the two systems become 
predominant for higher values of D, that is, for longer interaction times. This is evidently 
due to the cooperative interaction becoming more important for longer interaction times. 
Thus, we notice in Fig.l that < n > is in general higher compared to the one-atom 
micromaser for longer r. 

The variance in the cavity field is in general different in the two systems. It is generally 
very high near threshold, and hence, the sharp peaks in v at threshold appear at different 
values of D in the two systems as depicted in Fig. 2b. We also notice in Fig. 2b that 
the cavity field gets nonclassical properties {v < 1) for different ranges of D in the two 
systems. The photon distribution function also has different shapes as shown in Figs. 3, 
which clearly shows that the nature of the field depends clearly on the nature of the 
pump. For example, in Fig.3a, the field in the one-atom case is sub-Poissonian, while the 
field in the two-atom case is super-Poissonian. The field characteristics are vice-versa in 
Fig. 3c. However, the P„ in the present case do not indicate any existence of the so-called 
trapped states [4-6] . The Rabi frequencies in the two-atom Dicke system, represented by 
the eigen values Aq, A+ and A_ do not provide clear conditions for the trapped states as 
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one gets in the case of the one- atom micromaser [6]. There the conditions are essentially 
zeros of the function Sin-\/n + Ix, x being the dimensionless atom- field interaction time. 
The trapped states, are however, washed out by occasional presence of two atoms in 
the cavity of the one- atom maser [7,8]. Thus the present analysis gives an easy and 
transparent understanding of the disappearance of trapped states in multi-atom events. 

5. Comparison with one-atom two- photon micromaser 

The Dicke micromaser studied in the present paper has similarities with the dynamics 
of a one-atom two-photon micromaser [16] involving pumping of atoms, individually into 
a microwave cavity where the atoms make a two-photon transition [19] from the upper 
to lower level via an intermediate level. This is because the Dicke atomic system with 
J = 1 is equivalent to a one-atom three-level system with the middle level having equal 
frequency of separation uj from the upper and the lower levels. It may be noted here that 
uj is the transition frequency of the individual two-level atoms in the Dicke system. If 
the cavity eigenmode is in resonance with the two degenerate atomic transitions, and in 
addition, if they have dipole moments of equal strength, then the present theory can be 
employed to explain the dynamics of the two-photon cascade micromaser. In other words, 
the results in Figs. 1-3 can represent a one-atom two-photon cascade micromaser if we set 
the coupling constants of the two degenerate transitions gi — g\f2. The factor \f2 in the 
Dicke system represents its cooperative nature [14] . 

We need not stop at this comparison. A Dicke micromaser with pumping of three 
atoms at a time (all in their upper masing levels) will be equivalent to a degenerate 
five-photon cascade micromaser, and so on. 

However, a true two-photon process involves non-resonant one-photon transitions. If 
LO\ and ijJi are two atomic transitions, and v is the cavity mode frequency, then a true 
two-photon process should have large detunings l\\ — lo\ — v and A2 — v — loi^ that is. 
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Ai, A2 ^ gi- A two-photon resonance means Ai + A2 = 0, or c^i + c<j2 = 2z/. In the 
cascade two-photon micromaser discussed above, we have Ai = A2 = 0. An extensive 
comparison of cavity-QED of cascade two-photon with true two-photon processes can be 
found in [19] where it has been shown that the two dynamics are, in general, different. 

The present approach can also be applied to such a two-photon micromaser if we can 
accomodate the one-photon detuning in the derivation of the photon distribution function. 
This amounts to just setting the matrix elements = [Hi]^;?, = —A in Eq.4. For 

simplicity, we set gi — g\/2 here also. We present the photon statistics for the two photon 
micromaser for various values of A in Figs.4. We find that the operational characteristics 
are strongly dependent on A. We also notice that as A increases, the theory recovers the 
results for a two-photon micromaser derived by using an effective two-level Hamiltonian, 
a derivation and discussion of which can be found in [19]. 

For smaller values of A, we notice wriggles in the variation of < n > versus D in 
Fig.4a. These wriggles begin at a value of D[(x. r) and persist as it is further increased. 
The value of D at which the wriggles begin to appear increases with A. In other words, 
the smooth variations in < n > as in the case of an effective two-level system [17,18], occur 
for values of D having an upper bound which increases with A. Such characteristics have 
been discussed in detail in the context of the two-photon Jaynes-Cummings model in 
[19]. In addition to the above results, we also notice that the thresholds shift to higher 
values of D as A increases, and also, the range of D between two successive thresholds 
also increases with the one-photon detuning. The cavity field usually has sub-Poissonian 
photon statistics [v <1) between successive thresholds, and thus with a larger detuning, 
this model provides sub-Poissonian fields for wider range of D as shown in Fig.4b. 

6. Conclusion 

We have presented a theory for a two-atom Dicke micromaser and have brought out 
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the role played by the Dicke cooperativity on the micromaser action. Though at first 
glance it may look that the Dicke micromaser results mimic that for an equivalently 
pumped one-atom maser, the two systems are in general different. At places, however, 
the two systems have similar trends, an example of which can be seen in Fig.Sb. We 
further show that the micromaser dynamics involving the three-level Dicke atomic system 
is formally equivalent to a one-atom cascade two-photon micromaser. Interestingly, we 
note that the two one-photon coupling constants in the two-photon cascade transition 
has to scale -\/2 times of the one-atom coupling constants in the Dicke system for the 
two micromaser actions to have identical results. As mentioned above, this factor of -\/2 
originates from the cooperative nature of the Dicke system. With the introduction of a 
one-photon detuning in the two-photon process, the present approach also describes a 
one-atom two-photon micromaser. We have discussed in detail the role of the one-photon 
detuning in a two-photon micromaser action. Thus the Dicke superradiance effects can 
be quantitatively evaluated by scaling, as mentioned above, the results from the one-atom 
two- photon micromaser [16]. 

Thus the Dicke micromaser model discussed in this paper is not beyond the scope of 
experiments for verifying its results. The present paper gives the clue to verify results from 
micromaser actions involving pump mechanisms like the one studied here, which are rather 
difficult to achieve experimentally. This suggests a way to experimentally demonstrate 
the two-atom cavity-QED results [22] in a one-atom two- photon micromaser [16]. It may 
be recalled that the one-atom micromaser [1] demonstrated the phenomenon of quantum 
revival [23]. The same techniques may be employed in the case of the one-atom two- 
photon micromaser [16] and a proper scaling of its results should demonstrate the two- 
atom cavity-QED results [22] . 
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FIGURE-CAPTIONS 



Fig.l 

The cavity field intensity, proportional to < n >, versus the pump parameter D. 
nth = 0.1. N — 100 and 200 for the Dicke micromaser (a) and one-atom micromaser (b) 
respectively. This makes the total number of atoms that pass through the cavity in the 
two systems equal. The curve (b) is shifted upwards by 150 for clarity. 

Fig.2a 

Variation in < n > with respect to D for shorter interaction time. The parameters are 
same as in Fig.l. The full and broken curves are for the Dicke micromaser and one-atom 
micromaser respectively. 

Fig.2b 

Variation of v with D. The other parameters are same as in Fig. 2a. 
Fig.3a 

Photon distribution function for D — 2b. The other parameters are same as in Figs. 2. 
P{n) is sub-Poissonion {v — 0.60761) in the case of one-atom micromaser (broken curve), 
and is super- Poissonian {y = 1.27596) in the case of Dicke micromaser. 

Fig.3b 

P{n) for D — 50. The other parameters are same as in Fig.3a. For both the cases, 
P(n) is super-Poissonian at this value oi D. v — 1.09944 (broken curve) and v — 1.15871 
(full curve). 
Fig.Sc 

D — 400. P(n) is sub- Poissonian {v = 0.22911) in the Dicke micromaser system (full 
curve). One-atom micromaser (broken curve) provides a super-Poissonian {v = 1.05726) 
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in this case. 
Fig.4a 

< n > versus D for a two-photon micromaser for three different values of the one- 
photon detuning A. fith = 0.1 and N = 100. A = 100 (curve (a)), A = 150 (curve (b)), 
and A = 300 (curve (c)). For clarity, the curves (b) and (c) are shifted upwards by 40 
and 80 respectively. 

Fig.4b 

Variation of v with D for a two-photon micromaser. The other parameters are same 
as in Fig.4a. Here, the curves (b) and (c) are shifted by 4 and 8 respectively, for clarity. 
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